Effect of intrinsic instability of cantilevers on static mode Atomic Force Spectroscopy 
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We show that the static force spectroscopy curve is significantly modified due to presence of intrinsic 
cantilever instability. This instability acts in tandem with such instabilities like water bridge or 
molecular bond rupture and makes the static force spectroscopy curve (including "jump-off-contact" ) 
dependent on the step-size of the movement of sample stage. A model has been proposed to explain 
the data. This has been further validated by applying an electric field between tip and substrate 
which modifies the tip-substrate interaction. 
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Atomic Force Microscopy (AFM) has emerged as a 
powerful tool, having a wide variety of applications, from 
understanding the atomic level forces, Casimir force [![, 
friction at nanometer length scales @, H| to controlled 
manipulation of atoms jij. This relatively simple instru- 
ment has revolutionized our understanding of structures 
at the nanometer scale and hence the ability to manipu- 
late systems at atomic scales in a wide variety of subjects 
including material science, soft matter and biology [5, 6]. 

Since the time of its discovery by Binning et. al. Q, 
many attempts have been made to explain some of the 
non-intuitive features seen in these systems. For exam- 
ple, the force versus distance (f-h) curves § depend on 
whether the cantilever is approaching towards the sam- 
ple or retracting away from it (henceforth referred to as 
the " approach" and the "retract" , respectively) , leading 
to a hysteresis like behaviour as shown in fig03 The hys- 
teresis has traditionally been attributed to adhesion due 
to the layer of water existing on the surface of the sam- 
ple 0, [l(| , or rupture of molecular bonds 11 , [l2| , and 
has indeed been used to measure the "snap off" force. In 
this paper we show that the widely used practice of deter- 
mining the "snap off' force from the (f-h) curves can be 
erroneous because intrinsic instability in cantilevers can 
actually modify the (f-h) curves. We also show how one 
can properly interpret the (f-h) curves in the context of 
these instabilities. In an actual experiment the quantity 
measured is the cantilever deflection (d) as a function of 
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FIG. 1: (Colour Online) Schematic diagram of AFM tip and 
sample assembly (a) and deflection-distance curves at differ- 
ent step sizes taken on Si (b). The dotted line in (a) marks 
the equilibrium position of the cantilever in the absence of an 
external force, d is positive when measured upwards. The 
arrows in (b) show the direction of motion of the cantilever. 



the distance between the sample and the cantilever tip 
when the tip is in the equilibrium position (in the absence 
of any external force) (h). The force / = k c d, where k c is 
cantilever spring constant. It is important to note that 
although the steps in which the sample approaches or 
retracts from the cantilever (the z-controller resolution) 
is very small (« 0.025 A), the cantilever deflection d is 
only measured at discrete points in the whole path. In 
all our discussions below we will define "step size" (5h) 
as the distance between two such neighboring points, and 
assume that the distance between these points is covered 
smoothly without any noticable change to the deflection. 
If the maximum distance between the cantilever and the 
surface is h max and the number of data points acquired 
is A, then h max — N5h. In all our experiments, N has 
been kept fixed at 500 (in one direction) and hence Sh 
can be varied by varying h. 
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The two important parameters that one obtains from 
the experimental f-h curves are the "jump-into-contact" 
(JIC) distance obtained from the approach part and the 
"jump-off-contact" (JOC) from the retract part. The 
force, defined by /* = k c d* , where d* is the cantilever 
deflection at JOC, has traditionally been attributed to 
adhesion or molecular bond rupture. 

However, we observe experimentally that d* and h* 
(the tip-sample separation at JOC) depend on step size 
(Sh) as shown in fig. [2] In this paper we show both ex- 
perimentally and through theoretical analysis that these 
observed dependencies of d* and h* on Sh arise due to 
an intrinsic instability in the cantilever dynamics which 
manifests itself mainly due to the procedure of data ac- 
quisition in most AFM. We show, in particular, that both 
the instabilities (the intrinsic instability and the "snap 
off" instability) occur in tandem. We also find that in 
UHV-AFM where the "snap-off" instability is absent, one 
observes the JIC and JOC arising solely from the intrin- 
sic instability. The intrinsic instability arises mainly due 
to the motion of the cantilever in a non-linear force field 
and the two instabilities can be separated out in a real 
AFM experiment by acquiring data as a function of Sh. 
We support our inference by varying the tip-surface force 
and thus the instability in a controlled manner by apply- 
ing an electric field between them. We find that the shift 
of the observed d* and h* on Sh can be cleanly explained 
by our model. 

The experiments were carried out using a commer- 
cial AFM (Model CP II,Veeco) [13] using cantilevers 
(k c ps O.lN/m) with SiaN4 tip on cleaned Si wafers with 
natural oxide layer on it unless otherwise stated. The 
cantilever tip had a radius of curvature, Rt « 30 nm as 
determined by direct imaging. For experiment using ap- 
plied electric field, the substrate was gold film and a Si 
cantilever (fc c =0.2 N/m) with tip coated with Ptlr was 
used. A d.c. bias was applied to the tip from an exter- 
nal source and the sample was grounded. Experiments 
were carried out in a glove box with controlled Rh using 
flow of Ar gas at a temperature controlled environment 
at 28° C. The rate of data collection was 0.1 Hz for all 
the experiments presented here. Fig. [2] shows a set of h* 
and d* data plotted as a function of Sh. The data have 
been obtained from the typical (d-h) curves as shown in 
fig-ffl The data taken with two representative humidities 
(Rh=30% and 55%) are taken on Si surface with oxide 
(hydrophilic) . Another set of data (fig. [3]) are taken on a 
hydrophobic surface (created by etching the oxide layer 
using 50:1 (v/v) HF solution for 30 seconds). All the data 
show a definite trend. There are three regions in the data 
(barring the data taken on the hydrophobic surface). In 
region A, occurring at highest step size, we find that for 
Sh > Sh c i, both h* and d* , reach a limiting value which 
is independent of Sh. We call these limiting values h* min 
and d* nin , and they are almost independent of the Rh 
values. In the region C, that occurs for smaller Sh the 



—.22.0 

E 

C 20.0 
"° 18.0 
16.0 



d 



S 



-□- RH=55% 
-■- RH=30% 



1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 



28.0- 

— * 26 " 
£ 24.0- 
•S 22.0 - g 



: 20.0 



Step size (nrr 



4^ 



8h„. 



-O— RH=55% 
-•- RH=30% 



18.0- 
16.0 — 

1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 
Step size (nm) 

FIG. 2: Variation of d* and h* with step size. The closed 
squares and closed circles show the experimental curves for 
RH = 30% and the open squares and the open circles for RH 
= 55%. 



28.0 

26.0 

¥ 24.0 
c 

22.0 

TJ 

20.0 
18.0 
26.0 
24.0 

e" 22 

S 20.0 
*■= 18.0 
16.0 



-■-UHV-AFM data on Si 
4 — *— Data on Hydrophobic Si 

^^a^f^^. — •— Simulated data 



1.5 2.0 2.5 3.0 3.5 4.0 



Step size (nm) 

1 



1.5 2.0 2.5 3.0 3.5 4.0 

Step size (nm) 

FIG. 3: (Colour Online) Variation of d* and h* with step 
size. The closed squares show the experimental curves on Si 
taken in UHV-AFM, the closed triangles show the experimen- 
tal curves on hydrophobic Si and the closed circles show the 
simulated curves. The parameters used for the simulation are 
mentioned in the text. 



h* and d* again reach a limiting value h^ax an d d* max 
for Sh < Sh C 2 for hydrophilic surface. For Sh < Sh C 2 
both h* and d* become independent of Sh and h* max , 
d* max and Sh C 2 all depend strongly on Rh- In particular 
Sh C 2 is most sensitive to Rh and it increases as Rh is 
decreased along with the decrease in h* maxl d* max . For 
the hydrophobic surface (fig. [3}, there is no Sh c2 and h* 
and d* go on increasing as Sh is reduced. Data taken in 
an UHV-AFM is similar to that taken on a hydropho- 
bic surface (there is no 5h C 2). The data shown here are 
representative of a large number of data collected in the 
controlled experiment. In the region B which is the tran- 
sition region, both h* and d* increase as Sh is reduced. 
We propose that the two limiting regions in the data 
(region A and C) are due to the two instabilities that de- 
termine the cantilever motion. The instability at lower 
step size (region C) which depends on the humidity is 
due to the "snap-off" phenomena arising from the break- 
ing of the water bridge at the tip-substrate interface. A 
strong proof in favour of this is the observation that it is 
absent in the data taken on a hydrophobic surface and 
in the data taken using UHV-AFM. The instability at 
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higher step size (region A) is always present and arises 
due to the intrinsic instability that we describe below. In 
the following part we use a model proposed earlier by us 
14j to explain our observations. The motion of a can- 
tilever is modeled by a spring-ball system. The inherent 
nonlinearity of the cantilever due to its finite dimensions 
have not been introduced into our calculation, in order 
to keep things simple. Thus, we write the force balance 
equation for the static (or quasi-equilibrium) case as 



k c d = f ts {h + d), 



(1) 



The tip-sample interaction force ft s (h + d) is modeled 
by a combination of attractive van-der-Waals interaction, 
for a sphere plate geometry (which is close to the real 
situation) and the repulsive forces arising due to elastic 
interaction between the tip and the sample. In order to 
have a definite result Dejarguin-Muller-Toporov(DMT) 
force between the tip and the surface has been 
used. The tip-sample force is thus, formally, given by, 



fts(z) 




forz > a Q , 

\E*^R- t (a -z) 3 / 2 forz<a Q . 



(2) 



where H and Rt are the Hamaker constant and the ra- 
dius of curvature of the tip respectively. The attractive 
force is the only force present when h + d > the inter- 
molecular distance (ao), whereas when h + d < ao the 
force has a repulsive component, which inccrascs with 
reducing h. The repulsive component typically ensures 
that h + d > 0. It is interesting to note that, while the 
repulsive force is essential, the qualitative understanding 
of the f-h curves, comes even when the repulsive force is 
taken to be absent. In [14J we have ignored the repul- 
sive interaction for obtaining the exact solutions to the 
equation of motion of the cantilever. This will produce 
a slight deviation from the actual results, however, this 
will not change the conclusion. 

From eq. [1] and eq. [5] (in the region h + d > ao), we 
obtain the equation for the deflection (d) as, 



d(l + d) 2 + a = 0. 



(3) 



where d = d/h and a = H R t /6k c h 3 are dimensionless. 
The three exact solutions of this equation are already 
given in 0], therefore, we are not mentioning it here. 
For further discussion, we will refer these three solutions 
as di, d2 and d%. Solution d\ is real while d2 and ^3 
are either both real or complex conjugate of each other, 
depending on the parameters of the equation. Fig. 0] 
shows the simulated d-h curves (approach and retract) 
in presence of both attractive and repulsive part of the 
tip-sample interactions and also the analytical solutions 
of eq. [3] as a function of the tip-sample distance (h) for 
HRt=2.2 x 10~ 27 N.m 2 (appropriate for our experimen- 
tal conditions) and ao=0.172 nm. In simulated d-h curves 
we get realistic d\ (= d\.K) as in this case both attrac- 
tive and repulsive part of the tip-sample interactions are 
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FIG. 4: (Colour Online) Plot of the simulated d-h curves 
(black line) and the three analytical solutions (red symbols) 
of eq.[3]as a function of tip-sample distance (h) for the param- 
eters mentioned in the text. For analytical solutions, the open 
circles (di=di.h) and open triangles {dz=dz.h) represent sta- 
ble solutions and the open squares (d2=d2-h) represent the 
unstable solution. Three different regions (A,B and C) are 
shown in the inset. 



present. One can find out from fig.2]that d\ + h, obtained 
from simulation, is always positive whereas in case of an- 
alytical solution it is negative because the repulsive part 
of the tip-sample interaction has not been considered. It 
can be noted that the JIC position also matches quite 
well in these two cases. The solutions corresponding to 
di and d-$ approach each other and they meet at a point 
P (for example, at h ss 3.0 nm) and as h is reduced below 
this point they become complex. The distance of point 
P from the substrate (£) is the JIC point, which arises 
from the intrinsic cantilever instability 14]. Here di and 
d3 are stable solutions, while di is unstable. Hence, the 
tip will either equilibriate to the solution d\ or d3. 

We use fig. S] to explain the observed data. In most 
commercial AFM, during the process of data acquisition 
for the d-h curves, the motion of the cantilever is quasi- 
continuous, i.e., at each point the initial deflection (d) 
of the cantilever is determined by its final deflection at 
the previous point. In fig. [5] (in the inset), we show two 
examples of the paths traced by the cantilever (shown 
as steps). In one case (solid line steps), for relatively 
large step sizes (5h C 2 < Sh < 5h c \), the intrinsic instabil- 
ity dominates, and the jump from d\ to d3 occurs when 
the equilibrium position at the point just prior to the 
jump, takes the cantilever across di (marked by a cir- 
cle in fig. [4]). In the other case (dotted line steps), for 
relatively smaller step sizes (Sh ~ 8h C 2), the "snap off" 
instability dominates and causes a jump across the so- 
lution c?2. Here, lb rg determines the scale at which the 
water bridge snaps off, causing a jump across the solu- 
tion d2- If Sh > £, then, during the retract part, the 
cantilever tip will jump directly to the stable solution d3 
and d* and h* both become essentially independent of 
Sh. This corresponds to the region A, where the intrinsic 
instability is solely responsible for the JOC and which 
among other things depends on k c , Hamaker constant 
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FIG. 5: Variation of h* with step size as a function of bias 
voltage. The closed squares show the curve for V=0V, the 
open circles for V=2V and the closed triangles for V=4V. 
Variation of Sh c i as a function of applied voltage (V) is shown 
in the inset. 



H and tip radius R t . In region A we can thus identify 
£ w 5h c i. In region B both h* and d* increase as Sh 
is decreased. This is the region described above (solid 
line steps in fig. [4]). In absence of "snap off" instabil- 
ity, the region B extends all the way down to very small 
step sizes, as seen on the experiments on hydrophobic 
surfaces. On the other hand, if the "snap off" instability 
is present and Sh ~ 8h C 2, the JOC occurs when h ~ lf, rg , 
as discussed above (dotted line steps in fig. 2]). This is 
the region we identify as region C. In this region h* and 
d* are independent of Sh and lb rg is dependent only on 
Rh and Rt- Thus the qualitative discussion based on 
fig. [4] clearly identifies the regions of the observed curve 
and the instabilities that give rise to them. Thus in an 
actual experiment the (d-h) curves need to be taken as a 
function of Sh and the regions corresponding to the two 
instabilities can be clearly identified. In fig. [31 we show 
the data taken in an UHV-AFM, on a hydrophobic sur- 
face (which show only intrinsic cantilever instability) and 
the actual calculated curve based on the method above. 
The data for the curves taken with water bridges can also 
be fitted by taking h* nax as an adjustable parameter. It 
can be seen that the model proposed by us can explain 
the data in all the three regions. 

In this investigation, the "snap-off" occurs due to the 
instability of the water bridge that forms between the 
tip and the sample. 



It is shown earlier [17| that, for a 
sphere-plate geometry, depending on the radius of cur- 
vature of the tip, the water bridge configuration be- 
comes metastable for a particular sphere-sample sepa- 
ration when becomes ~ 1.0, where Rt is the radius 
of curvature of the tip. For the tip used i? t =30 nm, this 
should happen for h ~ 30 nm which matches very well 
with the value of h* nax ps 26.5 nm observed experimen- 
tally. The main proposal of the paper that there is an 
intrinsic instability of the cantilever can be further tested 
if we can modify the ft s in a controlled way. In region A, 
£ pa Sh c i and our model gives £ p=> ( 1,1 ^ Hiit ) 1 / 3 . From the 
experimentally determined H and R t we find that calcu- 
lated 5h c \ ~ 2.9 nm and experimentally obtained value 
is ps 3.0 nm. We have also checked that if we use softer 



cantilever having k c — 0.03 N/m (R t = 25 nm) then Sh c i 
shifts to 4.32 nm (calculated value of Sh c i in this case is 
4.29 nm). 

To establish the validity of our hypothesis we used elec- 
tric field to control ft s and obtained the (d-h) curves with 
applied electric field. The force due to the applied electric 
field [HI adds to the force term due to the van der Waals 
force and will give rise to an effective Hamaker constant 
H eff =H+C.V 2 {C w Constant). The enhanced H eff 
will make £ larger and will also shift Sh c i to a higher 
value. This can be clearly seen in fig. where we show 
the data taken in electric field. It has been observed that 
Sh c i clearly shifts to higher value with a small applied 
field. Since 6h cl w £ « C' 1 ™"* ) 1 ' 3 , in the applied elec- 
tric field one would observe that Sh c i ~ ~ V 2 / 3 . 
A plot of Shci 3 vs V 2 in the inset of fig. [5] shows that 
this dependence indeed exists. We also observe a shift of 
Sh C 2, h* max and d* max to higher values on application of 
electric field. This instability is related to water bridges 
[l^ | . It has been seen recently that water bridges become 
more stable in an applied electric field. The stability of 
the water bridge will be reflected in enhancement of Ibrg 
in an applied electric field leading to enhancement of the 
attractive force as shown in [l9j |. 

In summary, we have shown that the static d-h curves 
for an AFM, depends on the intrinsic instability of the 
microcantilever of the AFM. The phenomena like JIC 
and JOC occur even in absence of water bridge snap- 
off as in an UHV-AFM and on hydrophobic surface. At 
larger step sizes, the intrinsic instability dominates over 
the "snap off" instability, leading to erroneous results 
in the calculation of these forces. The instabilities due 
to "snap-off" forces dominate at smaller step sizes. We 
have also shown experimentally that the intrinsic insta- 
bility due to cantilever can be controlled by an applied 
electric field. 
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